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We propose a method to determine the topological structure of an event horizon far in the future 
of a spacetime from the geometrical information of its future null infinity. In the present article, 
we mainly consider spacetimes with two black holes. Although, in most of cases, the black holes 
coalesce and their event horizon is topologically a single sphere far in the future, there are several 
possibilities that the black holes do not coalesce eternally and such exact solutions. 

In our formulation, the geometrical structure of future null infinity is related to the topological 
structure of the upper end of the future null infinity through the Poincare-Hopf's theorem. Since 
the upper end of the future null infinity determines the event horizon far in the future under 
the conformal embedding, the topology of event horizon far in the future will be affected by the 
geometrical structure of the future null infinity. Our method is not only for the case of black 
hole coalescence. Also we can consider more than two black holes or a black hole with non-trivial 
topology. 



PACS numbers: 



I. INTRODUCTION 



We are often interested in the final state of a black hole 
spacetime after gravitational collapse. When matter col- 
lapses, it passes inside an event horizon and energy will 
be emitted to future null infinity in the form of gravita- 
tional waves. Since the amount of this energy is limited, 
we may expect that the spacetime will approach to a 
stationary state. If a black hole spacetime becomes sta- 
tionary far in the future, the final state of the black holes 
is considerably restricted. 

So far the final state of the black hole has been stud- 
ied by many authors and these studies are known as 
uniqueness theorem of black holes. Israel showed that 
the only static and topologically spherical black hole is 
Schwartzshild solution |Q or Reissner-Nortstrom solution 
||. In stationary axisymmetric situations, the unique- 
ness theorems of a vacuum black hole and a charged 
black hole were shown by Carter || and Robinson . [|| . 
Hawking proved that a stationary black hole have to be 
static or axisymmetric and that the spatial topology of 
the horizon is the only sphere |6) . 

In these works, however, the existence of endpoints of 
event horizon generators was not considered. In dynam- 
ical situation, the endpoints of event horizon generators, 
where the event horizon is indifferentiable, are very im- 
portant when we discuss the topology of an event hori- 
zon, because it causes the change of the spatial topology 
of an event horizon Q). This is revealed by one of au- 
thors |Q. When this was discussed in the topology 
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of an event horizon far in the future is assumed to be a 
sphere. On the other hand, the work of Chrusciel and 
Wald H implies that in a stationary spacetime, regard- 
less of the existence of a endpoint, the spatial topolo- 
gies of the connected components of black holes are only 
spheres under null energy condition. Therefore, if far in 
the future, the black hole spacetime settles down to a sta- 
tionary state, we expect the spatial topology of the event 
horizon will consist of some connected components with 
spherical topology. To realize the relation of these works, 
it is required to determine the topology of an event hori- 
zon far in the future. There is not the simple method to 
know topological structure of an event horizon far in the 
future, that is, to determine how many black holes the 
final state of the black hole spacetime will consist of. We 
will develop such a method in the present article. 

Since in an asymptotically flat spacetime JI^ ], (to 
be more exact, in a strongly asymptotically predictable 
spacetime | Q | ) an event horizon is defined as the bound- 
ary of the causal past of future null infinity, we expect 
the topology of an event horizon will be related to the 
topological structure of future null infinity. Because the 
topology of future null infinity in an asymptotically flat 
spacetime is S 2 x R, we may conclude far in the future the 
topology of an event horizon will always be a single S 2 . 
Nevertheless it is well known that there are some station- 
ary solutions including black holes which do not coalesce 
eternally. For example, they are Majumdar-Papapetrou 
solution ||, C-metric (l(J and so on. The spatial topol- 
ogy of the event horizon of such a solution is obviously 
not a single sphere. In these spacetimes, is the topology 
of the upper end of future null infinity really S 2 ? (Of 
course, since the future null infinity is open at the up- 
per end, this statement may be not well-defined.) If wc 
regard that the topology of the upper end of future null 
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infinity determines the topology of the event horizon far 
in the future, we may say that the topology of the upper 
end of future null infinity in such spacetimes is not one 
sphere but two spheres. 

In the present article, we consider that in such space- 
times with the black holes which do not merge eternally, 
the topological structure at the upper end of future null 
infinity pinches. Here we emphasize that its topologi- 
cal structure determines the spatial topology of an event 
horizon far in the future. In section [y, we will give the in- 
dex theorem and mention the relation between the topol- 
ogy of the upper end of future null infinity and geomet- 
rical structure of the null congruence of the generator 
on future null infinity. Since we discuss black hole space- 
times, we consider the only asymptotically flat spacetime. 
So we will give its definition in the rest of section [o|. In 
section we will discuss the method to examine the 
topology of the upper end of future null infinity. In sec- 
tion |y|, we will apply this method to C-metric as an 
example, to illustrate the detail concretely. 

Throughout this article, we use the abstract index no- 
tation as the component notation of tensors and it is 
denoted by Latin indices a, 6, • • • . 

II. PREPARATION 

How can we know the topology of an event horizon far 
in the future?? As mentioned above, an event horizon is 
defined as the boundary of the causal past of future null 
infinity, whose topology is S 2 x R 0|. So at first sight 
we think that the spatial topology of black holes is a sin- 
gle sphere finally. That is, all black holes coalesce and 
become a single far in the future. Is this real ? There are 
some exceptions to this. The examples of these excep- 
tions are Majumdar-Papapetrou solution and C-mctric. 
On the former, each black hole does not coalesce eternally 
by the force of electric repulsion. The latter describes two 
uniformly accelerated black holes connected by a string 
which are away from each other in the future. It is evi- 
dent that in these spacetime black holes do not coalesce 
far in the future. It seems that these are inconsistent 
with the topology of future null infinity, S 2 x R. 

Nevertheless, there is doubt that the topology of the 
upper end of ^ + is a sphere (FIG. 1(a)). For because the 
definition of is = J + (i°) — i , it is open to the 
future, where J + (i°) is the causal future of spatial infinity 
i°. Therefore, there is not evidence that the topology of 
the upper end of + is a sphere in such spacetimes. As 
shown in FIG. 1(b), we may expect that J'^ pinches at 
the upper end in the spacetime where two black holes 
does not eternally coalesce. 

First of all, we must mention the method to examine 
the topology of the upper end of future null infinity. The 
following corollary |ll| of the Poincare-Hopf 's theorem is 
useful. 

Theorem 1 Let M be a compact n- dimensional (n > 2 is 



(a) 




FIG. 1: The topological structure of the upper end of future 
null infinity in the case (a) where two black holes will coalesce 
in the future and in the case (b) where two black holes will 
not collide far in the future. In (a), ,J? + is a sphere at the 
upper end but in (b), J? + pinches at the upper end. 

an odd number) C r (r > 1) manifold with Si U E2 = dM 
and Si n S2 = 0. X is any C r ~ l vector field with at 
most a finite number of zeros, satisfying the following 
two conditions: (a) The zeros of X are contained in Int M. 
(b)X has inward directions at Si and outward directions 
at S2 . Then the sum of the indices of X at all its zeros 
is related to the Euler numbers 0/S1 and S2: 

X (E 2 )-x(£i) = 2index(A), (1) 

where index(A) is given by the alternating sum of the 
Morse number \ik as index = X)k( — l) k Mk- The Morse 
number ^ is the number of a critical point (zero of vector 
field X) whose index is k. The index of a critical point 
is given by the number of negative eigenvalue of Hesse 
matrix H ab = \7 a X b . 

This theorem means that when the topology of slices 
of a manifold changes, there must be a zero of the vec- 
tor field on it and equation (0) is satisfied. Because J? + 
is a three-dimensional null hypersurface, we can apply 
this theorem to Let us consider the tangent vec- 

tor field of the null generators on ^ + as the vector field 
in this theorem. Here, we should pay attention to the 
fact that, in this case, zero can exist only on the up- 
per end of J^ + , since is geodesically complete about 
own null generators by definition. Therefore we have to 
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be careful to apply the theorem. Nevertheless, it would 
be possible to relate the theorem of index and the con- 
gruence of null generators around zero. In FIG. 2 and 
the following, we temporarily consider small extension of 
future null infinity into its future directions for an expla- 
nation. From Morse's lemma, we suppose that the zero 
of the null generators is isolated. Since is three di- 
mension, the behavior of the congruence around zero can 
be classified and their index is determined. These are il- 
lustrated in FIG. 2. The left figures are three dimensional 
manifolds with index= ±1 zeros, and the right represent 
the behavior of the null congruence on the neighborhood 
around their zeros. From this figure, there will be type 
(a) zero on upper end of with index = +1. Without 
reaching to the upper end, the asymptotic behavior of 
shear and expansion tells us information about the zero. 
Then we expect shear dominates the congruence that is 
approaching to the type (a) zero, and expansion to the 
type (b)-(c) zero with index = — 1. If we assume that 
when two black holes do not coalesce, there should be a 
zero with index + 1 and we can say that the ratio of the 
shear to the expansion will diverge near the zero. 



a) Index=+1 




definition 1 A spacetime (M, g a b) is said to be asymp- 
totically flat at null infinity if there exists an unphysical 
spacetime (M, g a b) with g a b = 2 gab smooth everywhere 
except at the point i° where it is C >0 and conformal isom- 
etry \t : M — > ty(M) G M with conformal factor Q satis- 
fying the below conditions. 

1. J + (i°)UJ~(i°) = M-M. Thus, i° is spacelike re- 
lated to all points in M and the boundary of M con- 
sists off, J? + = j + (i°)-i° and J- = 

2. There exists an open neighborhood V of M = 
i U U such that the spacetime (V, g a b) is 
strongly causal. 

3. f2 can be extended to a function on all M which is 
C°° everywhere and C >0 at i°. 

4. On J?+ and J~ we have = and V Q ^ 0. 

5. The map of null direction at i° into the space of 
integral curves of n a = g ab V;,f2 on and J^^ is 
diffeomorphism. 

6. For smooth function on M — i with > on 
M U U which satisfies V a (uj 4 n a ) = on 

, the vector field Lo~ l n a is complete on J" + and 

J~. 

7. In the neighborhood of J? + and physical Ricci 
tensor behaviors as R a b = 0(fl 2 ). 

Let us note that there is gauge freedom in the choice of 
an unphysical spacetime [M , g a b) with an asymptotically 
flat physical spacetime (M, g a b) |L2]. This gauge freedom 
is most important in our discussion. 

If the spacetime (M,g a if) is an unphysical spacetime 
satisfying the above definition with conformal factor 0, 
so is [M, uj 2 g a b) with conformal factor 00 for the any 
function to which is smooth everywhere except at i° and 
positive everywhere. 



c) Index--1 





III. ANALYTIC METHOD 

In this section, we will investigate the relation of the ge- 
ometrical structure of future null infinity (the congruence 
of the null generators and curvatures) and the topology 
of event horizon far in the future. 



FIG. 2: The behaviour of vector field around the zeros of (a) 
index= +1 and (b),(c)index= —1. We can expect that in (a) 
the shear will become larger than the expansion and that in 
(b),(c) the shear will smaller than the expansion. 

Asymptotic flatness is essential in our discussion. We 
have to state the definition of asymptotically flat space- 
time before explaining the analytic method to examine 
the topological structure of J^ + . The definition men- 
tioned here is based on the reference |T3] . 



A. Physical Ricci Tensor and Unphysical Ricci 
Tensor 

The relation between the physical Ricci tensor and the 
unphysical Ricci tensor gives us a great deal of infor- 
mation. The physical Ricci tensor R a b is related to the 
conformally translated unphysical Ricci tensor R a b by 

Rab - Rab + 20- 1 V„V b 

+9ab9 cd (0- 1 V c V d - 30- 2 V c 0V d 0). (2) 
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By multiplying 2 and taking the limit — > we find 
that a vector n a = g ab \7 b must be extended smoothly 
to J? + and be null at because the first term vanishes 
from conditions 4 and 7 of the definition 1, the second 
vanishes from conditions 3 and 4, and the fourth vanishes 
from the condition 3,4 and g ab being smooth at J? + . 

n a is null at ^ + but is not null off of J? + and we put 
the shift from null vector as the following 

g ab n a n b = k (1) + k {2) 2 + ■■■ , (3) 

where kr{\ (i — 1,2, •••) are the functions independent 
of 0. The integral curves of n a on ,f + is null geodesic 
generators of but the integral curves of n a is not 
geodesies off of J r+ . Let us put the shift from the 
geodesic as 

n a V a n b -p {Q) n b = 0(0). (4) 

In particular, we put the contraction of the above equa- 
tion with l a = (d/d0) a as the following: 

l b n a V a n b = -p (0) + p ( i)0 + • • • . (5) 
From (|^) and (||), we obtain the relation between and 

1 3 

P(0) = 2 fc (!)' = " fc ( 2 )' P ( 2 ) = ~2 fc(3) ' ' " ' ^ 

where we used the relation of torsion free ; 

V a n h = V 6 n a . (7) 

B. Semi-Newman-Penrose Formalism 

Newman-Penrose formalism |l3| is the choice of the 
null basis which consists of a pair of real null vec- 
tors, l a ,n a and a pair of complex conjugate null vectors 
m a ,fh a . Now we consider n a = g ab \7 b as one of real 
null vectors in Newman-Penrose formalism. As we have 
mentioned before, n a is null on J? + but is not null away 
off J? + . Therefore, strictly speaking, this is not Newman- 
Penrose formalism. Nevertheless the Newmann Penrose 
equations are not affected by this difference. 

So the following orthogonality conditions are required 

l a m a = l a m a = n a m a = n a m a = (8) 
in addition to the conditions that vectors are null, 

l a l a = m a m a = m a m a = 0, (9) 
and the fact that n a = g ab V b is null only on J^ + , 

n a n a = k {1) + k {2) 2 + ■ ■ ■ . (10) 



We impose on the basis vectors the further normalization 
conditions, 

l a n a = -1, m a fh a = 1, (11) 

which is consistent with the definition of the vector n a 
and l a . Thus, the metric can be represented by 

9ab = (fc(i)0 + /c( 2 )0 2 H )l a l b -2l ia n b) +2m (a fh b) (12) 

in the basis vectors. Note that the first extra term of the 
above equation is caused by the equation ([l0|) . 

In order to investigate the behavior of the null 
geodesies on J? + , we have only to consider the values 
of spin connections and curvatures only in the neighbor- 
hood of J? + . Now we expand the each basis component 
of the equation ^ in the powers of conformal factor 0. 
The (n, I) component of the equation (^) gives 

O(0 2 ) = Qfc(i)-3p (O )+2 Ai(o) ^0- 1 

+ + 3p ( i) + 2fc (2) + 2/i (1) + k (1) qj + 0(0), 

(13) 

where we used the condition 7 of the above definition 1 
and the equations (H),(||) an d /U = ~m a fh b V a n b is the 
spin connection, which means the expansion of the inte- 
gral curves of n a . q is the function satisfying i°V a l b — 
ql b . 

So we obtain from the coefficients of the first and the 
second term 

5 

- 3p (0 ) + 2^( 0) = (14) 

and 

Rni + 3P(1) + 2fc (2) + 2 M(1) + k {1) q = 0. (15) 

Similarly, from other components we can obtain the 
relations between spin connections and Ricci tensors. 
From the (1,1) component, we obtain 

9-0, (16) 

4°' = 0. (17) 
From the (m, m) component, we obtain 

A(o) = A( ) = 0, (18) 

A(i) = R ( ±, (19) 

where A = —m a fh b V ' a n b is the spin connection which 
means the shear of the integral curves of vector field n a . 
The (n, m) component gives 

^(o) = V(o) = 0, (20) 

2v {1) - R(°l, (21) 
where v is defined by v = n a fh b 'W a n b . 



5 



C. Change of Conformal Factor 

Now it should be noted that there is gauge freedom in 
the conformal transformation considered above. Under 
the transformation — » 0' = 00, g ab — > <^ 6 = 2 g a b, 
we have 

n a -> n ,a = 0- 1 n a + &0~ 2 g ab V b 0. (22) 

In particular, on ^ + the above transformation becomes 

n a -> n' a = 0- 1 n a . (23) 

Under this transformation, the expansion of the null 
geodesic generators is transformed as 

^(o)-M , (o) = ^f "^ a V a (24) 

on Therefore, note that we can choose //m) as an y 

function independent of . Given fi', Q s , since the equa- 
tion ( |2~i| ) is merely an ordinary differential equation, there 
always exists satisfying the equation (£4j). From the 
equations do|)(|l4|), we find the relation 

M(o) = -P(o) = -\ k d)- ( 25 ) 

So from the equation (||), we see that in the gauge sat- 
isfying /i(Q) = 0, the null geodesies on J> + are affinely 
parameterized. In this article, we call this gauge affine 
gauge. When we discuss the upper end of we must 
not choose the affne gauge since in the affine gauge, there 
is not the upper end, that is, the null geodesies are com- 
plete in an affine parameterizing. Hence, we must choose 
the gauge such that ^ + is compact in the direction of n a , 
that is, at the upper end — > along the null geodesies 

on ,y + . 

From (|lS|), the zeroth order for the shear A, vanishes, 
which does not depend on the choice of the gauge; 

A' (0) = \o) = (26) 

holds. Under the transformation — ► 0' = 00, the first 
order is transformed as 

A 1 2 

A (1) - A' (1) = -f- - -rn*fh b V a V b - ^(m a V a0 ) 2 . 

From now on, ' represents the quantity after transforming 
from affine gauge into the gauge such that ,f + become 
compact. 

D. Null Generator Congruence and Weyl 
Curvature 

As mentioned before, the topology of the event hori- 
zon is expected to be related to the upper end of + . 



We discuss this in the context of asymptotic flatness. In 
an asymptotic flat spacetime, it is guaranteed that there 
exists a conformal embedding defined in the latter half of 
section D. The asymptotic flatness, however, accepts the 
further conformal transformation that is indicated by w 
as gauge freedom. Since the gauge transformation can 
become singular (uj can become zero or infinity) at the 
upper end of we should be careful to choose oj. The 
gauge freedom can makes the upper end of degener- 
ate into a point or take away to infinity in the direction of 
null generator tangent n a . Indeed, when we take a gauge 
choice in which n a becomes affine parameterized, the null 
generators are complete in affine parameterizing by defi- 
nition. Nevertheless, since we want to study the geomet- 
rical structure near the upper end of ^ + , it is important 
to choose a gauge 51' = ujfl in which the null generators 
are incomplete in affine parameterizing though is com- 
plete in a original parameterization n a — (d/du') a . This 
aspect results from the compactness of ^ + (or unphys- 
ical manifold) in the upper direction. Under this gauge 
choice, we discuss the topology of the upper end of . 

If we allow for uj to become angular depend irregularity 
at the upper end of , the irregular angular dependence 
also may change the topology of the upper end of ^ + . 
Then we do not accept such the dependence for the gauge 
transformation f2' = luQ from affine gauge. 

Newman-Penrose equationjl3| relates the shear on 
to Weyl curvature as the following. 

- n a V a \ + m a V a v = -(// + /2)A- (37-7)A 

+ (3a + + it - f)v + *( 4 )(28) 

We expand each term of the above equation in powers 
of and leave the only leading term. From ( |l8| ) the order 
of A is 

A = 0(0), (29) 

which is independent of the choice of gauge. The spin 
connection v can be written as 

v = -im a V a %)0 + O(0 2 ) (30) 

by the equations (||) and (||). We can choose satisfying 

m Q V Q %) = (31) 

from the relation (pa). For we can write m a as m a = 
m e (d/d9) a + im^(djd(j)) a using two spacelike vectors, 
(d/d9) a , (d/d(j)) a orthogonal to n a ,l a . Then if we trans- 
form uj from affine gauge into another gauge, from the 
equation (pi|) , the equation (||l]) becomes 

We may choose uj so that n a V a 0/0 2 does not depend on 
9, <j> and is positive everywhere. Therefore, we see that 
for satisfying (pl[), the order of v is 

v = O{0 2 ). (33) 
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Finally, let us compute the order of the spin connection 
7 defined by 7 = l/2(n a n b S7J b - n a rh b V a l b ). The real 
part of 7 is 

7 + 7 = n a n b VJb 

= /i (O )+O(0) (34) 

from the equations (^) and (|l4|). Since the imaginary 
part of 7 becomes 



7 - 7 = -n a m \7 a m b , 



(35) 



it depends on the direction of m a which we have not yet 
determined. Now let us determine the direction of m a 
as follows. We determine the direction of m a so that 
on ^ + , it will be parallelly transported along the null 
geodesic generators on That is, 



n a V n m b = on J?+. 



(36) 



Moreover in the direction away off it is parallelly 
transported along — l a , that is, 



l a \7 a m b = everywhere. 
Thus, 7 — 7 becomes 

7-7 = O(0). 



(37) 



(38) 



So we see that from the equations (|34|) and ( |38| ) , the order 
of 7 is 



7=^ + O(0). 



(39) 



From (p9|), ( p3[ ) and (p9[), we can write the equation 



n a V a X {1) = 3//( )A (1 ) - * 4 (i) 



(40) 



in the leading order. Integrating this, we obtain the re- 
lation 

A(i) =exp^3 J fi^du'j J du'9 iW C-3 J n^du' 

(41) 

where v! is the parameter of null geodesies on and 
is defined by n a = (d/du') a . If we transform conformal 
factor from affine gauge into the gauge such that J> + is 



compact and put = u' 



/X(0) is 



M(o) 



(42) 



So if the behavior of Weyl curvature on ,f + is ^4(1) ~ 
l/w' e as u' — ► 00, the ratio of the shear to the expansion 
is 



u l2 -'& + O(0' 2 ), 



(43) 



where e depends on u and is a function of a. Therefore, 
if there is a satisfying 2 — e > 0, this ratio can diverge, 



that is, the shear can be much larger than the expansion. 
This implies that pinches at the upper end. The 
point is that if Weyl curvature falls off later than some 
power of parameter of null generators, two black holes 
will not coalesce far in the future. 

Now that we obtain the sufficient preparation, we can 
restate our sufficient condition for not coalescing black 
holes as follows; 



Condition That Black Holes Do Not Coalesce 

When we transform the conformal factor from affine 
gauge into another gauge — > 0' = 00, the topology 
of future null infinity pinches at the upper end if there 
exists an unphysical spacetime (M, 2 g a b) with the 
conformal factor satisfying the below conditions. 

(1) There is a real number a such that = 0(9, (f>)u'~ a 
and lim u _KX3 — 0, where u and v! are parameters of 
the null geodesies on in affine gauge and in another 
gauge, respectively and 0(0. <f>) is the smooth function on 
a sphere which is positive and not singular everywhere. 

(2) There is a real number a > satisfying the following 
condition; For any number L > 0, there are positive 
numbers K and 5 such that if u' > K, < Q' < S, then 
\Q'u' 2 ~ e \ > L, where e is an exponent appearing in ^4(1) 
extended in the power of u' , that is, ^4(1) ~ u'~ e . 

Of course, in such an unphysical spacetime, the space- 
time is confomally embedded as illustrated in FIG. 1(b). 



IV. APPLICATION TO C-METRIC 

In this section we demonstrate the above condition by 
an example. Now we consider the vacuum C-metric as 
the example of the spacetime where two black holes does 
not eternally coalesce |Q , , [jl6| . The line element of 
the vacuum C-metric is written in the form of 

ds 2 = -Hdu 2 + H^dr 2 + 2Ar 2 H~ 1 drdx 

+ {G- 1 + A 2 r 2 H~ 1 )r 2 dx 2 + r 2 Gdtp 2 , (44) 

where 



H = 1 - 2mr 1 + 6Amx + ArG t . 
G = l-x 2 -2Amx 3 . 



A r G, (45) 
(46) 



BMS coordinate JL7| is expanded at large distance, but 
since in our discussion it is sufficient to consider the only 
neighborhood of we may transform the metric ( |44| ) 
to a BMS coordinate written by 



ds 2 



-fdu 2 - 2e 2f3 dudr 



-r 2 (e 



2 ^d0 2 



e 27 sin' 



2Udud9 

2a ^ 2 \ (47) 
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where four functions are 



2M{u,ff) 



f(u,r,9) = 1- 

U(u,r,6) = -(eg + 2c cot ( 
and M(u, 0) is given by 



c 2 («, 6 
4r 2 
c(u, 0) 



0(1 



o[ - 

r 



M u = -c 2 „ + -(c ee + 3c e cot0 - 2c), 



(48) 
(49) 
(50) 
(51) 

(52) 



We should note that these four functions are expressed 
by one function c(u, 9) whose differentiation by u is called 
the news function. We introduce the new coordinate by 



v — u + 2r. 



(53) 



In this coordinate, the metric ( [4T| ) becomes the following 
form. 

ds 2 = -(/ -e 2f3 )du 2 -e 2f3 dudv~2Udud8 

+-(v - u) 2 {e 2l d9 2 + e~ 27 sin 2 9d<p 2 ). (54) 

Furthermore, we introduce the new coordinate V = 1/v 
so that infinity along outgoing null geodesies will corre- 
spond to V = 0. The metric components in the new 
coordinate (u, V, 8, 0) are 



ds 2 



-;2B 



-(/ - e 2 ")du 2 + —dudV - 2Udud9 
+ ~ u) 2 (e 2 ^d6 2 + e~ 27 sin 2 d<i> 2 ). (55) 

In turn, we transform the above physical metric. 
dS 2 = -(/ -e 2/3 )0 2 du 2 + e 20 dud0-2U0 2 dud6 

+~(1 - 0u) 2 (e 2 ~<d9 2 + e- 2 ~< sin 2 ddcj) 2 ), (56) 

where the unphysical metric is related to the physical 
metric by g a b = 2 <lab and = 1/. 

Four functions (pq) , (fl9"l) , (|5fj|) and |[l]) are expanded as 

= -c 2 (u,6)0 2 + O(0 3 ), (57) 

7 = 2c(u,0)0 + O(0 2 ), (58) 

/ = l-4Af(u,e»)0 + O(0 2 ), (59) 
U = -(c<j + 2ccot0) + [-2c(c0 + 2ccot0) 

+2(2c „ + ice g + 4c 4 cot 9)]0 + O{0 2 ) (60) 

in powers of conformal factor 0. 

Now let us determine the basis vector mentioned in 



section [II from the metric section ([36|). It is easy to 
decompose the metric (|56| ) into the form of 

9ab = (fc(i)0 H )lJb - 2l (a n b ) + 2m(„m 6 ). (61) 



We see that the four dual vectors are 
n a = (d0) o , 



1 / 4C/0 2 

(1-U0) -e~ 7 ^.J du) 



2V2 V (i-u0) 2 

+e~<(d6) a -ie-~<sm6{d<f>) a V 



(62) 
(63) 



(64) 



(65) 



and we raise the indices of the dual vectors to obtain four 
basis vectors. 



1 / 4C/0 2 

(1-U0) -e" 7 - „, 2 W 



2\/2 v 'V (l-«0) 2 

+e 7 (d<9) Q -He^sin 



2e" 2 ^ 



o \ 

dv 



+ [4(/ - e 2f3 )e-^0 2 + lee- 4 * 3 - 27 



U 2 4 



(1-U0) 2, \d0 



d V 



-27-2/3 



u- 



2 



(l-u0) 2 V00, 
V2e~ 7 / .9 \ a . v^e 7 f d \ a 



r\ \ a 

1 90 



l-u0\d6j +l sm6(l-u0)\d(l) / 
_ a _ V2V 7 / d\ a , V2e 7 / g> \ a 
"' ~ l-«0 \ ~ Z sin(9(l - u 0) \d4>) ' 



(66) 
(67) 
(68) 
(69) 



From ( |62j ) and (66), we can compute the norm of n a as 
follows. 



»«»«„ = 4(/ - e 2f} )e-^0 2 + lee" 4 ^ 27 



[/ 2 4 



(1 -u0) 2 



(70) 



Since the first order in the equation ( f70| ) vanishes, we see 
that the choice of conformal factor = V is affine gauge. 
That is, in the choice of this gauge 



M(o) — -P(o) 



holds. The shear A in this gauge is 

A = — m a m v a nb 
= -4c u + O(0 2 ) 



3a 3 



m sm 



A 3 



-0 + O(0 2 ), 



(71) 

(72) 
(73) 



where we used the fact the news function of the vacuum 
C-metric Jlal is 



-3a m sin 8 



(74) 



4 A 3 u 4 ' 

and u is the affine parameter of the null generator of J? + . 
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We must transform conformal factor from affine gauge 
into the gauge such that ,J? + is compact and examine the 
ratio of the shear to the expansion of the null geodesic 
generators which plunge into the zero at the upper end of 
J^ + . In affine gauge, the shear in this direction behaves 
as u — ► oo and O — > 



A 



(75) 



in affine gauge. Under the transformation — > 0' = 00, 
there is the relation between affine parameter and the 
parameter of the geodesic in new gauge as follows 



Q+l 

log v! 



&(a = -1) 



where we set = 1/ v! a . 

We can obtain the ratio of the shear to the expansion 
from @. 



(77) 



where we used (pi|), (p7|), ( f75| ) and (|7q). So, we see that 
when u — > oc(u' — > oojTthis ratio diverges to infinity. We 
ensured the existence of the comformal factor 0' = 00 
such that the ratio of the shear to the expansion of the 
null generators on diverges. We also find easily 



lim = 0. 



(78) 



On the other hand, though Newman-Pcnrose formalism 
is ,of course, true, we confirm that the behavior of Weyl 
curvature is obtained directly from the metric 



1 



4(1) 



,ll-a ' 



(79) 



geometrical structure of the future null infinity. A suf- 
ficient condition for topologically non-trivial final black 
holes is proposed and applied to C-metric as an example. 
If some component of Weyl curvature has a late power 
law tail, since future null infinity pinches off at the upper 
end, two black holes will not coalesce far in the future. 
We make sure of the validity of our method on C-metric. 

While we concentrate on the gravitational radiation 
VP m of Weyl curvature as leading order contribution, 



there are possibilities that other component of curvature 
may contribute to the topological structure of future null 
infinity, for example, for Majumdar-Papapetrou solution. 
In Majumdar-Papapetrou spacetime, since gravitational 
(76) field radiation ^> u\ vanishes, we expect ^12) contributes 



to the topological structure of the upper end of in a 
higher order. 

In the present article, we mainly considered the 
case with several spherical black holes, since Wald and 
Chrusciel proved in a stationary spacetime the black 
holes are spherical under null energy condition. There- 
fore, if we expect the final state of a black hole spacetime 
become stationary, our method will give the number of 
black hole black holes with a spherical topology far in the 
future. 

Using the result of this article, what can we know 
about the collapsing stars from the observation of the 
gravitational radiation. The sufficient condition for topo- 
logically non-trivial black hole suggests the power of late 
time tail radiation will let us know about the final topol- 
ogy of black holes. Of course, to get rigorous statement 
about that, we should make clear the astrophysical mean- 
ing of the gauge choice used in this article. 



So if we choose uj so that 0<a<l, 0<e<2 holds. 
Therefore, two conditions (1) and (2) stated in the end of 
section |l| are satisfied. This implies that cross section 
of pinches at the upper end, as we have expected. 
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